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Abstract 

The effects of interpreting classical phase space distributions as 
Wigner functions, which is common in models of multiparticle pro- 
duction, are discussed. The temperature for the classical description 
is always higher than that for its Wigner function interpretation. A 
rough estimate shows that the corresponding correction is proportional 
to R~ 2 , where R is the radius of the interaction region, and that it 
is negligible for heavy ion scattering, but at the few percent level for 
e + e~ annihilations. 
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1 Introduction 



Much work is being done on the femtoscopy of the interaction regions. One 
of the main problems is to find the space-time distribution of the set of the 
freeze-out points, i.e. of the points where the hadrons are finally freed. This 
is known to depend on the momenta of the particles, which significantly 
complicates the problem. For reviews of the work in this field see e.g. [I], 

0,131. 

One can use several functions to describe the geometry of the interaction 
region in connection with the corresponding momentum distribution of the 
final state particles. The simplest is the classical phase space distribution for 
the particles at freeze-out F(p, x, t). Many models provide just that. This is 
immediately seen when classical equations are being used as e.g. the Euler 
equations from hydrodynamics or the classical Boltzmann kinetic equation. 
For a discussion of a number of cascade models from this point of view see [I] . 
The classical approach is intuitive and most useful to get a general picture 
of the situation. In principle it contradicts quantum mechanics, because it is 
not possible to ascribe to a particle simultaneously a position in space and 
a momentum. In practice, however, often the quantum corrections are not 
very significant. 

Another possibility is to use the Wigner function W(p, x, t). This is well 
defined in quantum mechanics. Its relations to the density matrices in the 
momentum and coordinate representations are 



W(K,X,t) = 
W(K,X,t) = 




where 



K=~(p + p'), q=p-p', X=i(x + x'), y = x-x'. (3) 

In a rigorously understood sense [5], [6] the Wigner function is the best 
quantum-mechanical replacement for the classical phase-space density. Heisen- 
berg's uncertainty principle is reflected by the inequality 
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|W( P ,X)| <7T 



-3 



(4) 



which follows from the definition of the Wigner function. Wigner functions 
integrated over momenta give the correct space distributions and integrated 
over the space give the correct momentum distributions. The quantum me- 
chanical averages of the type (x m p x ) cannot, in general, be reliably calcu- 
lated using the classical product with the Wigner function as weight, because 
they depend on the ordering of the noncommuting operators x and p x . The 
averages calculated with the Wigner function always give the quantum me- 
chanical average for the symmetrized (Weyl's ordering) product. E.g. 



Here and in the following the hats are used to distinguish operators from 
the corresponding classical quantities. The most annoying feature of the 
Wigner function is that only in very exceptional cases it is nonnegative. In 
fact, for pure states the Wigner function is nowhere negative if and only if 
the corresponding wave function is a Gaussian [7|. Fortunately, for mixed 
states this implies that any average over Gaussians satisfying (j3J) can be a 
Wigner function,which is enough to reproduce almost any shape, provided 
there are no peaks violating the bound (j3J). According to the class of mod- 
els described in the following section, in order to describe the multiparticle 
momentum distributions it is necessary to know the single particle density 
matrix pi(p; p'). As seen from (j2j), this can be calculated when the Wigner 
function is known. It cannot, however, be obtained directly from the clas- 
sical phase space distribution F(p,x, t). Therefore, models which yield the 
classical density usually tacitly assume that it is sufficiently similar to the 
corresponding Wigner function to replace it in formula (121). 

The purpose of the present paper is to study the relation between the func- 
tions F(p, x, t) and W(p, x, t). Our analysis suggests that this replacement 
is legitimate for heavy ion scattering, but overestimates the temperature of 
the system by several per cent for e + e~ annihilations. 

Still another possibility is to use the emission function [E], [H], [ID], related 
to the density matrix by the relation 



dxdp W(p, x)p 2 x x = -{p x x + 2p. 



» x xp x + xpl) 



(5) 
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where q and X are four-vectors and N is a constant factor. This formula 
is applicable for times after freeze-out has been completed. Then, in the 
interaction representation, the density matrix does not depend on time any 
more. The emission function is particularly convenient when the time spread 
of the freeze-out process is of interest, In the present paper only simultaneous 
freeze-out will be considered, so the emission function will not be needed. 

2 Simplifying assumptions 

The multiparticle system just after freeze-out is in some complicated, highly 
correlated state. Therefore, in order to deal with it, it is necessary to intro- 
duce approximations. The simplest would be to neglect all the correlations. 
Then the diagonal elements of an n-particle density matrix in the momen- 
tum representation, which is what one needs to get the n-particle momentum 
distribution, would be given by the formula 

n 

p nu (p 1 ,...,p n ;p 1 ,...,p n ) = JJpi( Pi ; Pi ), (7) 

i=i 

where u in the subscript stands for uncorrelated. In this approximation, 
however, for n identical mesons there are no Bose-Einstein correlations. Since 
the Bose-Einstein correlations yield important information about the particle 
distributions in coordinate space, a better approximation must be used. The 
next choice [TT] (for reviews see e.g. pQ, [2], [3]) is to introduce proper 
symmetrization over the momenta of identical particles. Then for n identical 
mesons 

n 

Pn(Pl J --- J P n ;Pl>---»Pn) = Cn^Ylp^p^Ppj). (8) 

P 3=1 

The summation is over all the permutations of the second arguments of p±. 
Symmetrizing also over the first arguments would just produce a constant 
factor n\, so there is no point in doing it. The normalization constant C n is 
now necessary to ensure the proper normalization of p n . With this choice, 
the single-particle and two-particle momentum distributions are 

P(p) = CVi( P ;p), (9) 
P(Pi,P 2 ) = C 2 (P(p 1 )P(p 2 ) + |p 1 ( Pl ;p 2 )| 2 ), (10) 
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where C n are normalizing constants and the hermiticity of the density matrix 
Pi has been used. Ansatz flS]) leaves out the final state interactions. For 
a study of e.g. resonance production this would be unacceptable, but for 
analyses of the interaction regions the available methods of removing the 
effects of final state interactions from the data are good enough [3] and Ansatz 
(JBj) is widely used. 

Let us assume further that, freeze-out for all the particles happens in- 
stantly and simultaneously at some time t = 0. With this assumption the 
emission function reduces to S(t)W(p, x) multiplied by a normalizing con- 
stant. Thus, there is no need to introduce an emission function besides the 
Wigner function, which greatly simplifies the discussion. Moreover, using the 
interaction representation one has, for t > 0, time independent density matri- 
ces and consequently time-independent Wigner functions. This assumption 
corresponds to a crude approximation. It would be better (cf. e.g. [3J) to 
assume that for each particle at its freeze-out its proper time r has some 
fixed value, common for all the particles. Then, however, the problem of 
comparison would become much more difficult. 

Finally, we assume that the particle density at freeze-out is given by 
the canonical distribution corresponding to some non-relativistic problem 
for noninteracting particles of mass m at temperature T in a force field 
corresponding to some potential V(x). The most characteristic implication of 
this assumption is that the space extension of the interaction region increases 
with increasing temperature. This is the case for most models, but it is 
not a law of nature. For instance, stars get hotter when they shrink. The 
canonical distribution is being used in many models. The assumption of 
a non-relativistic potential is not realistic, but it is sufficiently general to 
reproduce any size and shape of the interaction region, so it seems sufficiently 
flexible to provide qualitatively reliable results. 

3 The low- and high-temperature limits 

In the low temperature limit, classically, the particle rests (p = 0) at the 
minimum of the potential. Let us put there x = 0. Thus, 



Because of the sharp peak this cannot be interpreted as a Wigner function. 
In order to get a candidate Wigner function, F(p,x) must be smeared. 




(11) 
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The quantum-mechanical result is also easy to find. The particle must be 
in its ground state. Denoting the corresponding wave function ^ (x) we get 

W(p,x) = /(03^o(x+|)^o(x-|)e-^. (12) 

Now both position and momentum are spread around the point p = 0, x = 0. 
In the theory of fluctuations this effect is referred to as quantum fluctuations. 
Formula (jl2]) can be obtained by smearing ffTTT) . Therefore, smearing can be 
interpreted as a way of introducing quantum fluctuations. However, for each 
potential a different smearing prescription would be needed. Thus, at low 
temperatures the predictive power of the recipe: start with the classical 
distribution and smear it, is poor. 

According to our assumptions, the classical distribution is in general 

F(p,x) = iVe-^- /3V/ W. (13) 

Here and in the following iV denotes a normalization factor. Where the 
normalization factors are of no interest, we will use the same notation for all 
of them. The quantum- mechanical density operator is 

p = iV e -#-mx). (14) 

The difference between the classical and the quantum-mechanical expressions 
is that in the latter the kinetic energy does not commute with the potential 
energy. Let us note, however, that in the high-temperature limit [3 tends 
to zero. The commutator of the potential energy and kinetic energy terms 
in the exponent of (TH|) is proportional to [3 2 and, therefore, is negligible. 
Accordingly, in the high-temperature limit the two description are equivalent, 
as will be demonstrated more rigorously latter. 

The results from this section correspond to an effect which is well known 
from statistical physics. In the high-temperature limit the thermal fluctua- 
tions, common to the classical and quantum descriptions, usually dominate 
while in the low-temperature limit the quantum mechanical fluctuations, ab- 
sent in the classical case, are the important ones. 
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4 The smearing density operators 



When a classical phase-space distribution violates the bound PI, it is neces- 
sary to smear it. A way of doing this, is to introduce a (smearing) density 
operator as close as possible to the classical density distribution. Once a 
density operator is given, one can calculate from it a Wigner function which 
satisfies all the consistency conditions. 

Let us try, as smearing density operator, the operator 

p sm = Ne€ e -PVi*) e €. (15) 

The kinetic energy term has been split in order to make this operator her- 
mitian, as it should. The corresponding density matrix in the momentum 
representation is 

WintW) = iVe-^ (K2+ ^ 2) Jdx e-^*)-^, (16) 
and for the Wigner function one gets 

W sm (K,X) = Ne~^ Jdx e -mx) e -^(x-x )2 _ 

In the high temperature limit, /? — > 0, the second exponent in the integrand, 
taken with a suitable part of the normalizing factor, tends to 5(X — x) and, 
after integrating over x, one obtains the classical density distribution. For 
low temperatures, of course, no Wigner function can reproduce the classical 
distribution. 

Let us consider as example the harmonic oscillator with V^(x) = |fcx 2 . 
One gets 

W sm (p,x) = iVe-^-| fc * x2 , (18) 

where 

m* = m, k* = -k, (19) 

1 + rj z 

and 

V = \fa>- (20) 
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The parameter u = 




is the frequency of the oscillator. Note that with 



this smearing prescription the effective frequency u e ff = y depends on 
temperature. In order to get after smearing a distribution identical with 
F(p,x) one would have to make before smearing the substitution 

k -> (1 + r] 2 )k. (21) 

At high temperatures k* ~ k and the classical result is reproduced. At 
low temperatures, when 77 is large, k* <C /c and the x-distribution is smeared 
which avoids the contradiction to Heisenberg's uncertainty principle. 

As another example let us choose 

p sm = Ne -hm*) e -& e -hM*). (22) 
For the harmonic oscillator this yields again formula ffl8j) . but with 

m* = (1 + v 2 ) m, k* = k. (23) 

This time the smearing is in momentum space. A popular smearing prescrip- 
tion p], [BJ is 

psm = Jdpdx F(p,x,t)|V>(p,x))(^(p,x)|, (24) 

where the state vector \ip(p, x)) represents a bound state of one particle with 
positions close to x and momenta close to p. 

It is seen that various choices of p sm correspond to various smearing pre- 
scriptions. Each of them gives a reasonable Wigner function, but only with 
( I14p chosen as the smearing density operator the correct Wigner function 
is obtained. We will compare now in the general case the smearing density 
operator (fT5]) with the exact one (fT4"]l . 



5 Effective Hamiltonian 

The results obtained in the preceding section for the harmonic oscillator 
can be generalized to other potentials. One always finds that the smeared 
Wigner function corresponds to some Hamiltonian, but in general not to the 
true one for the system being studied. We will call this Hamiltonian effective 
Hamiltonian. It is defined by the relation 
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p m = Ne~ pA 'ff. (25) 
The smearing density operator (fT51) has the form 

Psm = Ne x e Y e x , (26) 
Thus, ignoring the irrelevant constants log N, 



-PH eff =\og(e x e Y e x y (27) 

The right-hand side can be evaluated from a simple extension of the famous 
Baker-Campbell-Hausdorff formula. The result is a series, in general infinite, 
of iterated commutators constructed from the operators X and Y . Since both 
X and Y are proportional to (3, this is a power series expansion in (3. An 
elegant and convenient method for calculating the coefficients of this series 

for a more general case, i.e. for log (e x e Y e^, has been described in ref. 
|15j . In our case an additional simplification occurs. Note that the operator 
e~ x e~ Y e~ x is the inverse of the operator e x e Y e x . Therefore, its logarithm 
equals +{3H e ff. On the other hand, the expansion of this logarithm can 
be obtained by taking the expansion for (|2Tj) and changing the signs of all 
the X-s and Y-s. These two prescription are consistent if and only if all 
the commutators with even numbers of factors have coefficients zero. For 
instance, for the smearing density operator (TT5]) one finds 



P 2 rr P 



H eff = H ~Y 



H, 



2 



+ •■■, (28) 



where H = + V(x) is the original Hamiltonian. The contribution of the 
single commutator vanishes as it should. 
In particular, for the harmonic oscillator 

Using a program in MATHEMATICA given in [15] it is easy to calculate more 
terms of this series. In order to get an effective Hamiltonian corresponding 
to the original phase-space density, and not to its smeared version, one must 
make in the Hamiltonian on the right-hand side the substitution (j2ip . This 
yields the Hamiltonian 
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H* = (l + -rfj H, (30) 

which reproduces, to second order in r\ the classical distribution ffTBl . 

For the harmonic oscillator, it is easy to compare directly, without using 
a smearing density operator, the Wigner function with the corresponding 
classical distribution. This is discussed in the following section. 



6 Classical density and Wigner function for 
the harmonic oscillator 

For a harmonic oscillator at temperature T, the Wigner function, or equiva- 
lently the density matrix, has been calculated by a variety of methods |16j . 
[6], [17]. The result is 

iy(p,x) = iVe-^ ( ^ + ^ x2) . (31) 
This is to be compared with the corresponding classical density 

F (p ;X ) = (^_y e -/W^+i^)_ (32) 

According to condition (jlj), if r] dass = \(3 c i ass uj > 1 the classicl density must 
be smeared before being interpreted as a Wigner function. The distributions 
(EU) and (|32]) coincide, if 

Vdass = tanh 7]. (33) 

At high temperatures, where r? and r] c i ass are both small, rj c i ass « r/ and there 
is no harm in interpreting the classical distribution as a Wigner function. At 
low temperatures, however, rj can be arbitrarily large, while rj c i ass never ex- 
ceeds unity. Then, interpreting the classical distribution as a Wigner function 
can lead to serious errors. 

An obvious question is, where on this scale are situated the temperatures 
in the range of some (100 — 200) MeV relevant for high energy scattering? 
The difficulty is that, what matters is the temperature in units of u, and u 
is not known. In order to get a rough estimate, let us make the admittedly 
crude assumption that the results for the harmonic oscillator can be used as 
a guide. For the harmonic oscillator 
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V 



This yields 



tanhr? = . — . (35) 

2 v /^) v /^fe) 



Choosing a value typical for high energy scattering, \/(7 2 (pi) = 300MeV, one 
gets as an approximation, which is very good when \t] — r/ c i ass \ is small, 

0.036 \ m 

1 — r T dass , (36) 



a 2 (xi) 

where <J 2 (xj) should be expressed in squared fermis. It is seen that for heavy 
ion high energy scattering, where typically ^ o 2 (xi) ~ 5fm, the correction 
is negligible. For e + e~ annihilations, however, where yj cr 2 (xj) can be as 
small as 0.7fm, the correction is about seven percent. We conclude that, 
interpreting the classical distribution as a Wigner function one always finds 
that the classical temperature is higher than the one corresponding to the 
Wigner function interpretation. Qualitatively this conclusion seems unavoid- 
able. The quantum fluctuations are not included in the classical description. 
In order to reproduce their effect it is necessary to increase the thermal 
fluctuations, which means increasing the temperature. The corresponding 
correction is probably negligible for heavy ion collisions, but may be at the 
few percent level for e + e~ annihilations. 

In order to obtain a Wigner function of the form ffTBl with f3 c i ass = (3, 
one has to start with the Hamiltonian 

V fr {■, , * 2 



H " = ^r n H ={ l+ 3" + -^ H ' < 37) 



which agrees with (130!) to second order in rj. 



7 Conclusions 

Numerous models provide classical phase space distributions for the parti- 
cles produced in multiparticle production processes. When describing Bose- 
Einstein correlations these densities, sometimes smeared, are being used as 
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if they were Wigner functions. Therefore, it is an interesting question how 
close, in situations encountered in particle physics, are the classical phase 
space distributions to their corresponding Wigner functions. 

Converting a classical phase space distribution to a Wigner function, when 
temperature is not very high, one should in principle consider quantum fluc- 
tuations. The simplest way is to assume that they are negligible. Our discus- 
sion, based on the analogy with the harmonic oscillator, suggests that this 
could be legitimate for high-energy heavy ion collisions, but probably not for 
e + e~ annihilations. 

In general, quantum fluctuations are negligible at high temperatures and 
important at low temperatures. For a given potential this means that, they 
are important when the interaction region is small, and unimportant when it 
is large. For the specific model discussed in the present paper, the correction 
goes like R~ 2 as seen from formula fl36l) . The correction always reduces the 
inferred temperature of the system. 

For potentials more complicated than that of the harmonic oscillator, it 
is convenient to perform the comparison of the classical phase-space distri- 
bution with the corresponding Wigner function in two steps. First one intro- 
duces a smearing density operator, which should provide a Wigner function 
easy to compare with the classical distribution. This is equivalent to the 
familiar smearing and yields a Wigner function which satisfies all the general 
consistency conditions. It can be done in an infinity of ways. Three are 
described in the text. The introduction of the smearing density operator is 
equivalent to the introduction of an effective hamiltonian which yields the 
same Wigner function as the smeared density operator. In the first two ex- 
amples discussed here, using the Baker-Campell-Hausdorff formula, one can 
obtain this effective Hamiltonian as a power series in (3. The leading term 
is the true Hamiltonian which confirms that in the high temperature limit 
(/3 — > 0) the crude estimate of the quantum fluctuations, as done by intro- 
ducing smearing, is good enough. This is implied by the fact that quantum 
fluctuations are negligible. 
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